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. The breakup of the 9 Be projectile on the 209 Bi target at bombarding 

energies above and near the Coulomb barrier is studied in the adiabatic two- 
center shell model approach. The effect of 9 Be— >n+2a breakup channel on 
complete fusion, elastic and inelastic cross sections is investigated. Results 
show that the breakup of the projectile 9 Be could be due to a molecular single- 



particle effect shortly before the colliding nuclei reach the Coulomb barrier. 
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)Q • I. INTRODUCTION 
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A. Review of the subject 

o _ 

The interaction of weakly bound or halo nuclei with stable targets at colliding energies 
around the Coulomb barrier is a very lively topic due to the increasing interest in nuclear 
reactions with radioactive beams. The existence and the role of the breakup process of 
weakly bound projectiles in complete fusion and scattering (elastic and inelastic) mechanisms 
have been extensively investigated in recent years both theoretically and experimentally 
|7Hl~7|, but there is not yet any definitive conclusion. 

There are contradictory theoretical works which predict either the suppression or 
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the enhancement of the complete fusion cross section due to the coupling of the relative 
motion of the colliding nuclei to the breakup channel. In a recent paper ||, coupled channels 
calculations for 11 Be+ 208 Pb have shown that the coupling of the relative motion to the 
breakup channel has two effects depending on the value of the bombarding energy, namely 
the reduction of the complete fusion cross sections at energies above the Coulomb barrier 
due to the loss of incident flux and the enhancement of the complete fusion cross sections 
at energies below the Coulomb barrier due to the dynamical renormalisation of the nucleus- 
nucleus potential. In |]-|3[], the effects of the breakup process on fusion were studied in the 
framework of both time-independent and time-dependent semi-classical coupled channels 
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theories using the imaginary part of a local dynamical polarisation potential associated with 
the breakup channel. In ||, the influence of the breakup process on fusion was studied 
by solving a three-body Schrodinger equation with a time-dependent wave packet method. 
The breakup process and its effects on fusion were described in in terms of inelastic 
excitations of the projectile to the continuum within the coupled channels formalism. 

In experiments of 9 ' 10 ' n Be+ 209 Bi J^To! and 9 Be+ 208 Pb |TJ, it was observed that the 
breakup process of both 9 Be and n Be projectiles suppresses the complete fusion above 
the Coulomb barrier. Reactions of 6,7 Li + 9 Be also reveal the same hindrance for fusion 
||, although contrary results also exist, e.g., 9 Be+ 9 Be |]17|| . This fusion hindrance was 
not observed in the reaction 9 Be+ 64 Zn ||16|| . We expect that the presence or absence of 
a suppression mechanism in fusion reactions involving 9 Be could be associated with the 
evolution of the last neutron of 9 Be in the approach phase of the collision before the nuclei 
reach the Coulomb barrier, which depends on the structure of the target or projectile. We 
also expect that molecular effects, which will be discussed later, do not occur in some 



reactions with Be. From an optical model analysis of the Be+ Bi elastic scattering [15 



the imaginary (absorptive) potential at the distance corresponding to the strong interaction 
radius increases with decreasing bombarding energy towards the Coulomb barrier energy, but 
this is not the case for the system 9 Be+ 64 Zn \TE\. This behaviour of the imaginary potential 



indicates that strong absorption channels are still open for bombarding energies very close 
to the Coulomb barrier and it was associated with the Coulomb breakup of 9 Be. Moreover, 
it was pointed out that other absorption channels related to orientation and deformation 
effects of 9 Be cannot be excluded. Such effects have been reported for the scattering of 9 Be 
by 40 < 44 Ca and 39 K 



B. Object of the work 

The aim of this paper is to study the breakup of the 9 Be projectile on the 209 Bi target 
at energies above and near the Coulomb barrier in the adiabatic two-center shell model 
(TCSM) jl9[ approach, as well as its role in the complete fusion and scattering processes. 
The TCSM is based on the assumption that the nucleons can be described with molecular 
states during the heavy ion reaction. This means that the motion of the nuclear centers has 
to be adiabatically slow compared with the rearrangement of the mean field of nucleons. 
Whether molecular orbitals are formed depends on the ratio of two characteristic times, 
the collision time r c ~ 2R/v = 2R(2E / M A)~ 1 ^ 2 and the nuclear period or single-particle 
relaxation time r s ~ 2R/vFermi = 2R{2e F jM)~ 1 l 2 , where E is the laboratory bombarding 
energy, A the projectile mass number, and e F the Fermi energy in the target or projectile 
nucleus with valence nucleons. According to the ratio of r c / r s ~ 3-4 for the studied reaction, 
we expect that molecular orbitals and the polarisation of valence nucleons by the field of 
the other nucleus may have time to develop. 

We study the neutron level diagram in the approach phase of this collision in order to see 
possible single-particle effects (nucleon transfer or single-particle excitations). Since the 9 Be 
has a small separation energy (S n = 1.67 MeV), the breakup channel ( 9 Be-^n+ 8 Be) could 
be related to a diabatic single-particle motion of the last neutron of 9 Be into the continuum 
|f20|| or to a neutron transfer to 209 Bi. Since the nucleus 8 Be has no bound states, as the 
Coulomb field of the target ( 209 Bi or 210 Bi) in its vicinity is very strong and the relative 
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motion is considered adiabatic, we assume that the decay 8 Be— > 2a occurs automatically 
and immediately. This assumption is experimentally supported by the observation of a 
particles originating from the decay of 8 Be |[21|| . The capture of a part of the projectile after 
the breakup is associated with the incomplete fusion channel which will not be considered 
in the present work for the sake of simplicity 

In general, the neutron promotion from 9 Be to 209 Bi or into the continuum could be 
related to three types of transitions [22"] between the molecular adiabatic single-particle 
levels: the transition induced by the radial relative motion of nuclei (radial transition ~ 
d/dr), the transition induced by the rotation of the internuclear axis around a fixed axis 
in the laboratory system (rotational transition ~ jl , j and / being the single-particle 
total angular momentum and the total angular momentum of the system, respectively), 
or the transition induced by the change of the orientation of the intrinsic symmetry axes 
of deformed nuclei, which is dependent on the angular momenta of the nuclei (~ I[, I' 2 
with I[ + I' 2 — I — L and L being the orbital angular momentum of the collision). In the 



TCSM [19] used in the present work, only radial and rotational transitions between the 



molecular adiabatic single-particle levels can occur because the intrinsic symmetry axis of 
nuclear shapes lies parallel to the internuclear axis. For arbitrary orientations of the intrinsic 
symmetry axes of deformed nuclei, a more general TCSM |23| would have to be used. 
The radial transition occurs between adiabatic states with the same angular momentum 
projection j z and is largest at the point of an avoided crossing (pseudo-crossing), while the 
rotational transition occurs between adiabatic states <fi a , (fib with j% — j\ ± 1 and is largest at 
the crossing point of levels. Moreover, the radial excitation mechanism describes the non- 
adiabaticity of the relative motion (Landau-Zener mechanism ||24|| ) and grows with increasing 
radial velocity, the diabatic single-particle motion being a limit for large radial transition 
probabilities ]|20| . In the diabatic limit of the single-particle motion [2D|, the nucleons do 



not occupy the lowest free single-particle levels as in the adiabatic case, but remain in the 
diabatic levels keeping their quantum numbers during a collective motion of the nuclear 
system. 

The study of single-particle level diagrams |?2 can be useful to understand some reaction 
channels and to simplify coupled channels calculations. A microscopic formulation of the 
excitation and nucleon transfer in the TCSM and its role in the scattering process of nuclei 
with loosely bound nucleons was presented in [2j| on the basis of the molecular particle-core 
model. A similar formulation was suggested in [2£| in the framework of the orthogonalised 
coupled-reaction-channel theory. Molecular single-particle effects are experimentally well 
established phenomena |22| in heavy ion collisions and have been reflected, e.g., in large 



structures observed in elastic and inelastic excitation functions of reactions between light 
nuclei. 

In sect. 2, trajectories, complete fusion, elastic and inelastic cross sections are calcu- 
lated in the framework of a classical model, using the potential energy obtained with the 
adiabatic TCSM along with phenomenological friction coefficients. The role of the breakup 
of 9 Be on complete fusion, elastic and inelastic cross sections is incorporated by using the 
Landau-Zener approach for the breakup induced by the radial mechanism. We compare the 
calculated results with experimental data in sect. 3. A summary and conclusions are given 
in sect. 4. 
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II. MODEL 



A. Adiabatic TCSM 

The basic microscopical model for a dinuclear system is the TCSM. The TCSM describes 



nuclei, namely the united nucleus and the two incoming nuclei. In the TCSM |19j used in 
the present work, the nuclear shapes are defined by a set of coordinates: 

1. The relative distance r between the centers of mass of the two nuclei (or the elongation 
A of the system, or the distance (z\ — z 2 ) between the centers of the two oscillators). 

2. The neck parameter e = Eq/E' is defined by the ratio of the actual barrier height Eq 
to the barrier height E' of the two-center oscillator. The neck grows with decreasing e. 

3. The transition of the nucleons through the neck is described by the mass (charge) 
asymmetry r\ = (A 2 — Ai)/A (rjz = [Z 2 — Z^/Z), where Ai(Zi) and A 2 (Z 2 ) are the 
mass (charge) numbers on both sides of the neck and A = A\ + A 2 [Z = Z\ + Z 2 ). 

4. The deformations /3j = di/bi (i — 1, 2) of axial symmetric fragments are related to the 
ratio of their semiaxes. For separated nuclei, the semiaxes aj and 6j can be related to 
the parameter of quadrupole deformation by using the known expansion of the nuclear 
surface in spherical harmonics. 

In the adiabatic approach, the parameters of the momentum-independent part of the 
single-particle TCSM potential are determined by the assumption of conservation of the 
volume enclosed by the nuclear surface, continuity of the potential and its derivatives be- 
tween the two centres and by the requirement that the TCSM potential barrier with height 
E is located at z = 0. The parameters related to the momentum- dependent part of the 
potential (spin-orbit and I 2 - like terms) are determined by some interpolation between the 
appropriate values for the compound nucleus and for the separated fragments for large r. In 
the present work, the same interpolation method as in [19] is adopted, but the parameters 
k and fi of the Nilsson model for the spin-orbit interaction are used |27j . The nuclear-radius 



constant is r = 1.2249 fm and the oscillator quanta are hum = 41 • A~ 1/3 MeV. 

The single-particle levels, characterised by the projection j z of the total angular mo- 
mentum, are given in the harmonic oscillator potentials only up to constants V associated 
with the depths of the two potential wells. In order to compensate this disadvantage, we 
use three parameters related to the depths of potential wells of the compound nucleus Vo 
and the separated nuclei V°° for large r. These parameters are expressed by the sum of 



experimental nucleon separation energies [28— 30| and corresponding Fermi levels. For finite 



relative distances, we interpolate the values of the depth of potential wells V, as 
v= { Vo + (Vr - V )^ if R neck > 1 

1 1 VP°, if Rneck =0 J ' 

where w % z are the frequencies of the oscillator along the internuclear z axis for the configu- 
ration of the system with a radius of the neck R nec k = 0. The u)\ and u° z are the z-oscillator 
frequencies for the nuclei in the united system and for the compound nucleus, respectively. 
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It is necessary to point out that are added to the single-particle energies taking care to 
identify every level in association with its localisation in only one of the two potential wells 
for separated nuclei. 

In order to avoid unrealistic shapes for small relative distances r (compact shapes) and 
a large mass asymmetry 77 of the colliding nuclei in this version of the TCSM, it is assumed 
that the mass asymmetry r\ decreases, e.g. linearly |19|], with decreasing relative distance 
r from the touching configuration of the nuclei. An adequate shape parametrisation for 
such highly asymmetric compact configurations was used to study cluster emission in the 
superasymmetric TCSM [31]. Since the motion of the system is mainly considered in the 



entrance channel up to relative distances near the touching configuration, the parametrisa- 
tion used for compact shapes does not influence the results. Moreover, compact shapes of 
the system do not occur in the approach phase of a collision for bombarding energies above 
the Coulomb barrier owing to the repulsive core of the nucleus-nucleus potential caused by 
the centrifugal effect and by diabatic P2| , |35| or structural forbiddenness effects ]3JJ . 

The orbital angular momentum independent part of the total potential of the two-center 
system is calculated with the macroscopic-microscopic method of Strutinsky as 

V = E macr (r, e, 77, r)z, (3i,fo) + 5E sheU (r, e, rj, r) Z , fa, (3 2 ) + 5E pair (r, e, r], r] Z , f3 u (3 2 ) - V C n- (2) 

The smooth macroscopic part of the potential energy E macr is expressed by the sum of 
the Coulomb energy E c and the nuclear surface energy E N as in the liquid drop model, 
and they contain a part related to the binding energy of the system and another one related 
to the nucleus-nucleus interaction potential. In order to take into consideration the finite 
range of the nuclear force and the diffuse nuclear surface, we use the nuclear surface energy 



En obtained by the Yukawa-plus-exponential method |3J| with shapes generated by the 
TCSM potential. The fluctuating microscopic shell and pairing corrections, 8E sheU and 
5E pair , are obtained using the single-particle levels of the TCSM and the Strutinsky method. 
The dependence of these microscopic corrections on temperature is neglected because the 
dissipated energy for the studied reaction is only a few MeV |T5j . The value of the potential 
(0) is normalised to the potential energy Vcn of the spherical compound nucleus. 



B. Classical trajectories 

Classical trajectories for a system with reduced mass /i, total energy E and orbital 
angular momentum L = {xr 2 d(p/dt in the center of mass system, are obtained by solving the 
generalised Lagrange equations in the relative distance r and orientation = {r, tp}) of 
the two centers of the nuclei 



d 
It 



d£ 



d(dqi/dt) 



d£ 

— = Fi(q,dq/dt) = ~2_,lij f j {q)dq j /dt, 



(3) 



where £ 
model 



36] 



T — V is the Lagrangian and F{ are the friction forces. As in the surface friction 
we assume a diagonal friction tensor 7^ = c®5ij (radial and tangential friction) 



with friction form factors fj(q)- The form factor of the friction forces is calculated as 



f TtV = [V(E N - E N )f 



(4) 
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Ejsf being the nuclear surface energy for the configuration of the system with a radius of the 
neck R nec k = 0. We have assumed in (Q) that frictional effects start with the overlapping 
of the two nuclei. Theoretical justifications for the classical form of the friction forces 
were given in []36[| , where the perturbation theory was used to study the coupling between 
the collective and the intrinsic degrees of freedom for small overlap between the nuclei. 
The tangential part of the friction is considerably weaker than the radial part [|36|| . The 
phenomenological friction forces in the classical trajectory express the irreversible coupling 
between the relative motion of nuclei and the degrees of freedom that are not explicitly 
considered (e.g., collective vibrations of 209 Bi), for small overlap between the nuclei outside 
the touching distance. Moreover, this phenomenological friction does not take into account 
the reversible coupling (diabatic single-particle motion due to the Landau-Zener mechanism) 
between the relative motion and the intrinsic degrees of freedom [S7|, which is much stronger 
for large overlap between the nuclei inside the touching distance. The friction parameters 
c° r and c° are usually considered as free parameters and their values are determined from 
extensive analyses of fusion cross sections and inelastic scattering data for many systems and 
bombarding energies throughout the periodic table p6| , p8| . For simplicity, we assume that 



the reduced mass fx and the friction coefficients c° do not depend on the temperature of the 
system. Moreover, the slight increase of the reduced mass with decreasing relative distances 
for weakly necking shapes around the touching configuration is neglected f39fl . In f40fl , it 



was pointed out that the inclusion of this contribution does not lead to any appreciable 
change of the observables in the exit channel. For large overlap between the nuclei, these 
transport coefficients can be calculated in an adiabatic approach within the linear response 
theory jyj. In the present work, classical trajectories are calculated using the radial and 
the tangential friction coefficients c° = 4 x 10~ 23 MeV -1 s and c° = 0.01 x 10~ 23 MeV -1 s, 
respectively p6 . 



C. Complete fusion, elastic and inelastic cross sections 

The solution of the classical equations of motion ([|) leads to captured and scattered (elas- 
tic and inelastic) trajectories. Since the initial mass asymmetry rj of the system ( 9 Be+ 209 Bi) 
is larger than its Businaro-Gallone mass asymmetry [Q, we assume that the complete fu- 
sion channel is mainly related to the relative distance r between the centers of the nuclei. 
In general, the complete fusion cross section can be expressed as 

a fus = °"ca P (E, I) ■ Fcw(E, I), (5) 

where a cap , Pqn and l cr are the partial capture cross section, a hindrance factor for complete 
fusion and the highest orbital angular momentum (a sharp angular momentum cut-off model) 
for trajectories which reach the critical distance r cr , respectively. At the critical distance 



r cr , the shells of the individual nuclei dissolve into those of the compound nucleus p? 
Expression (||) is valid if the mechanism hindering the complete fusion occurs before or after 
the Coulomb barrier has been passed. The partial capture cross section a cap is defined as 

^ = ^(2Z + 1).T,(E), (6) 
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where E is the bombarding energy. The transmission coefficient 7] through the Coulomb 
barrier Vbi with orbital angular momentum L = Ih is approximated by that of a parabolic 
barrier with angular momentum dependent frequency hui 

1 

In (|5p, the factor Pqn takes into account other hindrances for the compound nucleus 
formation. In the reaction 9 Be+ 209 Bi, we assume that the hindrance is mainly related to 
the breakup of the projectile 9 Be (e.g., 9 Be— m+2a): 

P CN = l-Pbu P (V,l), (8) 

where Pb up is the breakup probability. The hindrance factor (H) was used in and 
suggested experimentally in 0,0]. The role of Pb up in the complete fusion cross section 
depends on the existence of the breakup channel which is derived, in the present model, from 
the diagram of neutron levels calculated with the TCSM. Since for orbital angular momenta 
< / < l cr the breakup probability Pb up changes much less than the partial capture cross 
section a cap , we use an average value Pb up with respect to orbital angular momenta. The 
complete fusion cross section (W) is factorised as 



nh 2 ^ 21 + 1 

2^E ^ 1 + exp[2n(V Bl - E)/hu l } 



;i-n„ P (E)). (9) 



For simplicity, we introduce the approximations [[EJ hu\ = hcu and tbi = tb, with huo 
and tb being the frequency and the position of the I = Coulomb barrier Vb, respectively. 
Moreover, the summation in Eq. (|9]) is replaced by an integration and the term in brackets 
yields an analytical expression as the formula for the complete fusion cross section in the 



Glas-Mosel model [43]. The Glas-Mosel model is used, instead of a dynamical fusion model as 



the dynamical energy surplus model |44j], to calculate complete fusion cross sections without 
breakup because we found that the radial kinetic energy loss of nuclei ( 9 Be+ 209 Bi) before 
reaching the Coulomb barrier is negligible compared with the static Coulomb barrier V B , e.g., 
< 0.1 MeV at the bombarding energy 1.15Vb (E c . m . = 46 MeV) for < I < l cr . However, 
frictional effects on the radial motion of nuclei are important in calculating the breakup 
probability Pb up because this probability exponentially depends on the radial velocity of 
nuclei as we will discuss further on. 

In the complete fusion process, the probability for a breakup process induced by the 
Landau-Zener mechanism at an isolated pseudo-crossing of adiabatic single-particle levels 
can be expressed as 

Pbu P = Plz-, (10) 

and in the scattering process as 

P bup = P% + (I - P%) ■ Ptz\ (11) 

since the system passes through the pseudo-crossing point twice, first entering (e) and then 
leaving (Iv ) the interaction region. Here, we have made a distinction between (e) and (Iv ) 
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because the radial velocity changes due to frictional effects. The first term in Eq.flTTD is 
Plz, instead of P^ z ■ (1 — P^), because of our assumption of an immediate breakup of 8 Be 
— > 2a after the neutron in 9 Be has been removed from the nucleus entering the interaction 
region (the neutron does not get bound again with 8 Be when the nuclei leave the interaction 
region). At an isolated pseudo-crossing of adiabatic single-particle levels Ex, E 2 for a relative 
distance r , the Landau-Zener probability P3J for a nucleon transition between the single- 
particle levels is 

P L z = e-^ G \ (12) 

with 

G° = z n , } H ' 12 ' 2 , , , (13) 



where | H' 12 | is a coupling matrix element between the diabatic states ex, e 2 and is equal to 
half the closest distance between the adiabatic levels E\,E 2 at r . Diabatic states €i,e 2 are 
a set of states which run smoothly through the crossing. The radial velocity v°, which is the 
only quantity concerning relative motion in the above expression, can be calculated from the 
classical equations of motion (|3|). Since the radial velocity v° depends on the incident energy 
E and the orbital angular momentum I, P^z depends also on E and I. The Landau-Zener 
formula is derived under the following assumptions: H[ 2 is smaller than the relative kinetic 
energy of the nuclei, (ei — e 2 ) is a linear function of r, and v° is constant. The restriction 
to constant velocities implies that the point of avoided level crossing is far from the turning 
point of the relative motion of the two nuclei. The assumption of a constant relative velocity 
can be approximately fulfilled in a region where the nucleus-nucleus potential is shallow. In 
general, the neutron transition at a pseudo-crossing can be calculated semiclassically f22| 



using a time-dependent Schrodinger equation for the neutron wave function expanded in 
a basis of two diabatic wave functions related to the diabatic levels ex,e 2 and a classical 
equation for the radial relative motion of the nuclei. 

The scattering function 9(1) is calculated from the classical trajectories. Elastic and 
inelastic scattering cross sections a are obtained from the classical scattering cross section 
a d as 



da_ _ b(9) 
dQ sin 9 



db{9) 



d9 



■ P sc , (14) 



where the classical scattering cross section a c \ is given in terms of the scattering angle 9 and 
of the impact parameter b = 1%/ ^2fiE. The hindrance factor for scattering P sc = 1 — Pb up 
(Pbup is obtained from Eq. (|TTD) takes into account the fact that the number of scattered 
projectiles is reduced by the breakup process. We define as elastic events those non-captured 
trajectories which show less than 0.5 MeV energy loss, and as inelastic ones all other non- 
captured trajectories. 



III. RESULTS OF THE CALCULATION AND DISCUSSION 

Since the collision above the Coulomb barrier is discussed, we first consider 9 Be and 209 Bi 
as spherical nuclei (ratio of the semiaxes (3i = 1), and then study the deformation effects. 



S 



The assumption of spherical nuclei is suitable if the change of the nucleus-nucleus potential, 
due to orientation and deformation effects of the nuclei, is small in comparison with the 
surplus of bombarding energy above the potential. It is known p8[ that the deformation 



effects in the entrance channel may be important for collisions very close to the Coulomb 
barrier, where they mostly cause an increase of the complete fusion cross section. 

In order to obtain a nuclear shape for the touching configuration of the nuclei similar 
to the one supplied by the overlap of the two nuclear frozen densities [32|J3,J2|, the neck 
parameter e should be set at about 0.75. With this value of e, the neck radius and the 
distance between the centers of the nuclei are approximately equal to the corresponding 
quantities in the dinuclear system formed by the overlap of the two frozen densities. 



A. Neutron level diagram of the adiabatic TCSM 

Fig. 1 shows the neutron level diagram of the TCSM for the system 9 Be + 209 Bi — > 218 Fr 
as a function of the relative distance r between the spherical nuclei. For r pa 5.5 fm (the 
relative distance between the centers of mass of the two hemispheres of the spherical com- 
pound nucleus), we recognise the neutron single-particle states of the spherical compound 
nucleus 218 Fr and for large r the neutron single-particle states of the spherical separated 
nuclei 9 Be and 209 Bi. The quantum numbers l,j for these asymptotic neutron levels are 
shown. The single-particle levels quickly change in the vicinity of r ~ 5.5 fm because of 
the abrupt change of the nuclear shape from a spherical shape (at r pa 5.5 fm) to compact 
asymmetric shapes with a large value of the neck parameter e = 0.75 (for r > 5.5 fm). Small 
values of the neck parameter e (0 — 0.2) are more realistic for such compact shapes near 
the spherical shape because the potential energy of the system decreases with decreasing 
values of e |[22|| . The nuclei 9 Be and 209 Bi only reach spherical shapes for large r due to 
a polarisation effect caused by the smoothing of the barrier which joins the two oscillator 
potential wells ||19| . This evolution to spherical shapes depends slightly on the value of the 



neck parameter e. In Fig. 1, L* and H* denote the highest occupied level for the nuclei 9 Be 
and 209 Bi, respectively. The large gap for the neutron subsystem of 209 Bi associated with the 
neutron magic number N=126 is observed. The internal arrow (A) indicates the distance 
r t = 9.82 fm corresponding to the touching configuration of the nuclei and the external one 
(C) indicates the relative distance r = 13.15 fm for the configuration of the system with a 
radius of the neck R nec k = 0. The Coulomb barrier (B) is located between these two arrows 
at tb = 11.4 fm. 

In Fig. 1 at r ~ 12 fm, we can see two very close pseudo- crossings (denoted by 1 and 
2) between the highest occupied state of 9 Be with j z = 3/2 and two unoccupied states of 
209 Bi with j z = 3/2. Pseudo-crossings (denoted by 3 and 4) between these states also occur 
at relative distances close to the touching configuration of the nuclei. The level denoted by 
L corresponds to a level of 9 Be with j z = 1/2, occupied by two neutrons. This level shows 
a pseudo-crossing (denoted by 5) at r pa 12.5 fm with an unoccupied state of 209 Bi with 
j z = 1/2 as well as a crossing (denoted by 6) with another state of 209 Bi with j z = 3/2. In 
general, at pseudo-crossings 1-5 transitions induced by a radial mechanism can occur. The 
rotational mechanism can induce a transition at crossing 6, although this type of transition 



is much weaker than the radial one 22 . The radial transitions 1-4 are associated with the 



most favourable breakup channel Be — >n+2a (Q = —1.57 MeV). 
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In the present work, we will only consider transitions at the pseudo-crossings 1-2 of Fig. 1. 
The last neutron of 9 Be could be promoted into the continuum at pseudo-crossing 1 and could 
be transfered to 209 Bi at pseudo-crossing 2. At the pseudo-crossings 1-2, the applicability 
of the Landau-Zener approach is quite good. The values of | H[ 2 | and 4-(ci — e 2 ) at the 
pseudo-crossings 1-2 are 0.1073 MeV, 0.0184 MeVfm" 1 and 0.1250 MeV, 0^0352 MeVfm^ 1 , 
respectively. It could be expected that transition probabilities at pseudo-crossings 3-4 are 
much smaller than at pseudo-crossings 1-2, because the radial velocity decreases at the 
touching configuration of the nuclei due to strong frictional effects. 

From the neutron level diagram, one can observe that the shells of individual nuclei 
dissolve into those of the compound nucleus between 7.5 — 10 fm, which agree with the 
empirical critical distance r cr = 8.01 fm obtained by Galin's formula [f5J. We will use the 



value r cr obtained by Galin's formula as the critical distance. It is important to note that 
the crossing and the pseudo-crossing of levels occur nearly at the same relative distances, 
independent of large variations in the parameter sets of the TCSM. 



B. Nucleus-nucleus potential 

Fig. 2a) shows the total nucleus-nucleus potential for zero orbital angular momentum 
of colliding nuclei as a function of the relative distance r between the spherical nuclei (solid 
curve). We compare the results to those obtained using a total nucleus- nucleus potential 
(dashed curve) which is calculated as the sum of experimental binding energies j|7| of the 



nuclei (Bi + B 2 — Bi 2 ) and a sudden nucleus-nucleus interaction potential. The B\, B 2 
and B\2 are the experimental binding energies of the colliding nuclei and the compound 
nucleus, respectively. The sudden nucleus-nucleus interaction potential is obtained with the 
double-folding method using the Skyrme-type effective density- dependent nucleon-nucleon 
interaction and a realistic two-parameter symmetrised Woods-Saxon function for the density 



of nuclei [[42]]. Microscopic corrections to the double- folding potential are included in the 
binding energies. 

For large relative distances and for relative distances around the barrier, the potential 
calculated with the adiabatic TCSM is similar to the potential obtained with the double- 
folding method and small differences are explained by different nuclear shapes used in the 
two approaches. The height and the frequency Tiu of the barrier for both total nucleus- 
nucleus potentials are about 26 MeV and 2 MeV, respectively. For relative distances smaller 
than the position of the barrier, the adiabatic TCSM-potential decreases slower than 
the double-folding potential because of a slower decrease of the macroscopic nuclear surface 
energy (Fig. 2b). Fig. 2a) also shows the adiabatic TCSM-potential around the barrier for 
the prolate deformed nuclei in the pole-to-pole configuration (dashed-dotted curve), where 



the static deformations (ratio of the semiaxes) of Be [45[ and Bi are taken as f3\ = 1.2 



and @2 — 1-05, respectively. The deformations also take the effects of polarisation of the 
nuclei into account, induced by the mean-field of the partner-nucleus. In this case, the value 
of the barrier decreases by ~ 1.3 MeV, with approximately the same frequency hu, and its 
position is shifted to tb ~ 11-5 fm. The role of this potential for complete fusion very close 
to the barrier will be discussed further on. 

In Fig. 2b), the nuclear surface energy is normalised to the nuclear surface energy 
for the configuration of the system with a radius of the neck R nec k = 0. The small increase 
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of this energy at r fa 12.5 fm is related to the formation of a neck between the colliding 
nuclei. The double-folding potential shows a repulsive core for small relative distances which 
reflects the action of the Pauli principle in the Skyrme-type effective density- dependent 
nucleon-nucleon interaction. The numerical values of this sudden potential are similar to 
the diabatic potential even though both potentials are conceptually and physically not 
equivalent. The repulsive character of the diabatic potential is mainly related to diabatic 
particle-hole excitations. We would like to stress that there is practically no difference 
between the adiabatic and diabatic potentials of the two nuclei for relative distances larger 
than the one corresponding to the touching configuration, because the number of pseudo- 
crossings of single-particle levels is very small [|32lp3 . 



C. Complete fusion cross section and breakup process 

Fig. 3 shows the complete fusion cross section (jf us as a function of the inverse of the 
bombarding energy in the center of mass system 1/E c m .. Experimental fusion cross sections 
(full dots) for 9 Be + 209 Bi are from [ITJ. It was pointed out in ]TUj that the experimental 



fusion data for 9 Be + 209 Bi include the complete fusion of 8,9 Be. In ||10|| , it was assumed 
that the 8 Be dissociation after the breakup of 9 Be (— > n + 8 Be) occurs in a time much 
longer than the reaction one. However, we assume the opposite that the decay 8 Be— > 2a 
occurs automatically and immediately once 9 Be breaks up before reaching the Coulomb 
barrier. The capture either of 2a particles by the target 210 Bi (in the case of the neutron 
transfer) or of all projectile fragments (n + 2a) by the target 209 Bi (in the case of the 
neutron single-particle motion into the continuum) can also occur and, therefore, contribute 
to experimental complete fusion cross sections. We do not include these events because this 
is in itself a complex problem that needs to be investigated further. The ocurrence of these 
events as well as the incomplete fusion process depend on both kinematical aspects of the 
reaction (e.g., the angle between the 2a particles, which depends on bombarding energy) 
and dynamical aspects (e.g., motion of 2a particles in the mean field of the target). Hence, 
the present calculations yield a lower limit for complete fusion cross sections. Comparing 
this lower limit with the experimental complete fusion cross sections, we expect that the 
capture of all individual components of 9 Be by the target, after the 9 Be breakup, is not 
very significant. Our results are also compared with the experimental complete fusion cross 
sections (full triangles) obtained for the similar reaction 9 Be + 208 Pb [12]. 



In Fig. 3, we can see that the theoretical complete fusion cross sections calculated with 
expression (^]) (dotted curve) agree well with the experimental complete fusion cross sections. 
The adiabatic TCSM-potential for the spherical nuclei is used for these calculations and the 
results weakly depend on the value of the critical distance r cr . Predictions (without breakup) 
of the Glas-Mosel model for this nucleus-nucleus potential are shown by a solid curve. The 
theoretical hindrance factor (1 — -Pfe up (E cm .)) for complete fusion changes with increasing 
bombarding energy from 0.76 (24% of breakup) for bombarding energies 1.05 times the 
value of the Coulomb barrier to 0.66 (34% of breakup) for bombarding energies about 1.15 
times the value of the Coulomb barrier. These values agree with the experimental value of 
0.75 obtained for the studied reaction [ TO | and with the experimental value of 0.68 ± 0.07 



reported in |12| for the similar reaction 9 Be + 208 Pb. The calculation of complete fusion cross 
sections with breakup (dotted curve in Fig. 3) stops at the bombarding energy 1.05 times 



11 



the value of the Coulomb barrier for two reasons, namely the assumption of spherical nuclei 
is not suitable for bombarding energies very close to the Coulomb barrier and the breakup 
probability (|9|) cannot be calculated with the Landau-Zener formula ( |T2"D for subcoulomb 
trajectories because a radial velocity of the nuclei is not defined in the classically forbidden 
region. The arrow in Fig. 3 indicates the value of the Coulomb barrier of the adiabatic 
TCSM-potential for the spherical nuclei which is about 40 MeV. The value of the Coulomb 
barrier is obtained by subtracting (B\ + B2 — -B12) from the value of the barrier of the total 
nucleus- nucleus potential of Fig. 2a). 

In Fig. 3, the dashed-dotted curve shows the complete fusion cross sections calculated 
with the Glas-Mosel model (without breakup) for the adiabatic TCSM-potential with prolate 
deformed nuclei in the pole-to-pole configuration. The complete fusion cross section for 
bombarding energies very close to the Coulomb barrier can increase due to static deformation 
and orientation effects which lower the value of the Coulomb barrier to ~ 38.7 MeV. Due 
to the deformation effects and to the orientation of the intrinsic symmetry axes of deformed 
nuclei, the height of the Coulomb barrier can be larger or smaller than the one obtained 
with spherical nuclei. The complete fusion cross section for bombarding energies very close 
to the Coulomb barrier can be considered as an averaged value for different orientations of 
the intrinsic symmetry axes of the deformed nuclei. It could be expected that the role of 
the breakup (Fig. 4, discussed further on), suppressing the complete fusion ([|), is reduced 
for decreasing bombarding energies towards energies below the Coulomb barrier, and an 
enhancement of the complete fusion can be a result of large transmission coefficients for 
the smaller fusion barriers associated with certain orientations (e.g., pole-to-pole) of the 
deformed nuclei. 

Fig. 4 shows the breakup probability Pf, up as a function of the orbital angular momentum 
I of captured trajectories (0 < / < l cr ) for two values of the bombarding energy E c m , 
namely 42 MeV (dashed curve) and 46 MeV (solid curve). The breakup probability Pb up 
is calculated as the product of two independent Landau-Zener transitions (|T^) at the very 
close pseudo-crossings 1-2, respectively, of Fig. 1. For < I < l cr and E c . m .= 42 MeV 
(46 MeV), the breakup probability P bup decreases by a factor 1.3 (1.7) whereas the partial 
capture cross section o cap (|6|) increases by a factor 6.5 (10.5). For a fixed bombarding 
energy E c m , the breakup probability decreases with increasing orbital angular momentum I 
because the potential energy increases and the radial velocity (kinetic energy) decreases at 
the pseudo-crossing points. For a fixed orbital angular momentum I, the radial velocity at the 
pseudo-crossing points and, therefore, the breakup probability increases with an increasing 
bombarding energy. Results are similar to those obtained when the double-folding potential 
determines the trajectories. 

D. Elastic and inelastic scattering cross sections 

For a bombarding energy E c m . = 46 MeV, Fig. 5 shows the scattering function 8(1) 
(upper part), the radial distance r min of closest approach of the nuclei (middle part), and 
the dissipated energy Edi SS (lower part) as a function of the orbital angular momentum 
I for non-captured trajectories, using two different potentials. For small orbital angular 
momenta I, the TCSM results (solid curve) are quite similar to those obtained using the 
double-folding potential (dot-dashed curve). Differences mainly appear for large values of / 
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due to the different tails of the nucleus-nucleus potential. For the adiabatic TCSM-potential 
with spherical nuclei, we can observe that inelastic events are associated with low orbital 
angular momenta (12 — 22h) corresponding to angles smaller than the rainbow angle 9 r 
p2]| (where the scattering function 9(1) has a local maximum). The dissipated energies for 
these trajectories are a few MeV. The results do not depend strongly on the value of the 
radial friction coefficient c , (Fig. 6, upper part). Since most of these inelastic trajectories 
show a turning point inside the radius of the pseudo-crossings 1-2 of Fig. 1, their classical 
inelastic cross sections are reduced by the effect of the breakup of 9 Be. For trajectories with 
I = 12 — 16h, the classical turning point is relatively far from the position of the pseudo- 
crossings 1-2 (r « 12 fm) and we use the Landau-Zener approach to calculate the breakup 
probabilities Pb up - In Eq. (|TT|), the transition probabilities P^\ an d P^z are calculated as 
in the complete fusion channel. 

Fig. 6 (middle part) shows the angular distribution of inelastic events. The classical 
inelastic angular distributions depend weakly on the radial friction coefficient c°. They 
slightly change due to small shifts of the rainbow angle 9 r between 90 — 95° corresponding 
to dissipated energies between 1.5 — 3 MeV for the scattered trajectories with the rain- 
bow angle 9 r . Classical inelastic cross sections are very large in the neighbourhood of the 
rainbow angle (rainbow scattering) and the cross section becomes infinite at 9 r . Since the 
projectile 9 Be has no excited bound states, the inelastic events are related to the excita- 
tion of the target 209 Bi. The breakup of 9 Be reduces the classical inelastic cross sections 
by 29 — 34% for scattering angles smaller than the rainbow angle 9 r , but their values still 
remain appreciable. The experimental point is extracted from |15| and corresponds to the 
maximum of the angular distribution of the inelastic collective (vibrational) multiplet 3~ 
of 209 Bi at around 2.6 MeV for a scattering angle of 95°. This multiplet could not be 
observed at angles smaller than ~ 90°, which was explained by an increasing background 
most probably originating from the tail of the elastic peak. The large values of the elastic 
cross section (one-two orders of magnitude larger than the inelastic one) for angles smaller 
than ~ 90° support this claim (Fig. 6, lower part). This reason could also explain that all 
the other 209 Bi inelastic channels were observed with negligible cross sections [[15] . 



We have assumed that a corresponding classical trajectory is associated with the peak 
of the experimental angular distribution for this inelastic channel. Moreover, the classi- 
cal trajectory corresponds to the mean trajectory. In fact, the width of the experimental 
angular distribution for a given inelastic channel is caused by quantal effects and statisti- 
cal fluctuations |[38| . Although the experimental point of Fig. 6 (middle part) cannot be 
reproduced in the present classical treatment, we would like to note that the classical tra- 
jectory corresponding to the dissipated energy of 2.6 MeV has a scattering angle very close 
to the experimental one being in the vicinity of the rainbow angle 9 r . In reality the classical 
singularity at the rainbow angle is removed due to quantal effects, statistical fluctuations 
and also experimental broadening [ 3l| . For example, the inelastic angular distribution in 
the semiclassical approximation for rainbow scattering |22] using the Airy formula shows a 
smooth exponential decrease from ~ 200 mb to ~ 100 mb with increasing scattering angles 
around the rainbow angle 9 r . However, the semiclassical method of scattering [22 is not 
appropriate for the reaction studied here because the reduced mass fi of the system is very 
small and the value of de Broglie wavelength h/y / 2fiE of nuclei in the relative motion is 
comparable to their nuclear radii. The experimental value in Fig. 6 (middle part) could be 
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on the "dark side" (8 > 9 r ) of the angular pattern of rainbow scattering. This experimental 
value of the inelastic cross section can, however, well be reproduced by the coupled channels 
method ||15|| , where the relative motion of the nuclei is treated quantum mechanically. 

Fig. 6 (lower part) shows the angular distribution of elastic events, normalised to the 
Rutherford cross section or. For smaller scattering angles, the elastic cross section obtained 
with the adiabatic-TCSM potential slightly deviates from the Rutherford cross section due to 
the deviation of the tail of the adiabatic TCSM-potential from the Coulomb tail for large r. 
The increasing values around ~ 90° occur because the scattering angles are near the rainbow 
angle 8 r . The effect of the breakup in the elastic scattering is not considered because the 
turning point for elastic trajectories is outside the radius of the pseudo-crossings 1-2 of Fig. 
1. Results obtained when the double-folding potential determines the elastic trajectories are 
closer to the experimental data. 

For decreasing bombarding energies, the reduction of the projectile penetration leads 
to a decrease of the dissipated energy for inelastic events. We found that the rainbow 
angle 8 r in the scattering function 8(1) disappears for E c m m 44 MeV. The scattering angle 
for inelastic trajectories (E^ss > 0.5 MeV) becomes larger than 95° increasing to larger 
angles with decreasing bombarding energies. Inelastic cross sections become smaller with 
decreasing bombarding energies because the increase of dd(l)/dl is larger than the increase 
of the prefactor of (|14]) proportional to E"^ . Moreover, the breakup probability Pb up of 9 Be 
decreases with decreasing bombarding energies for inelastic trajectories (E^ ss > 0.5 MeV) 
where the radial distance r m i n of closest approach is smaller than ~ 12 fm (inelastic events 
with smaller values of /) . This probability changes with decreasing bombarding energy from 
0.34— 0.29 (/ = 12 — 16K) for bombarding energies 1.15 times the value of the Coulomb 
barrier to 0.2 — 0.11 (/ = 8 — 10h) for bombarding energies about 1.05 times the value of 
the Coulomb barrier. As in the complete fusion channel, the effect of the breakup of 9 Be on 
inelastic cross sections is diminished with decreasing bombarding energies. 

IV. SUMMARY AND CONCLUSIONS 

The breakup of the 9 Be projectile on the 209 Bi target at bombarding energies above and 
near the Coulomb barrier has been studied in the adiabatic two-center shell model approach 
as well as its effect on the complete fusion and scattering processes. The neutron level 
diagram of the two-center shell model reveals two very close pseudo-crossings between the 
state j z = 3/2 of the last neutron of 9 Be and two unoccupied states of 209 Bi with the same 
projection of the single-particle total angular momentum j z , shortly before the colliding 
nuclei reach the Coulomb barrier. A radial velocity-dependent transition of the last neutron 
of 9 Be into the continuum and to 209 Bi at these pseudo-crossings of molecular levels leads 
to the most favourable breakup channel of 9 Be ( 9 Be — > n + 8 Be). Since the nucleus 8 Be 
has no bound states, the Coulomb field of the target ( 209 Bi or 210 Bi) in its vicinity is very 
strong and the relative motion of the nuclei is considered adiabatic, we assumed that the 
decay 8 Be — > 2a occurs automatically and immediately. Breakup probabilities depending on 
the bombarding energy and the orbital angular momentum have been calculated with the 
Landau-Zener approach and trajectories obtained from a classical model using the adiabatic 
TCSM-potential along with phenomenological friction coefficients. 

The effect of this breakup channel on complete fusion and scattering (elastic and inelastic) 
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cross sections was studied. From the value obtained with the Glas-Mosel model without 
breakup, the complete fusion cross section reduces by 34% for the bombarding energy 1.15 
times the value of the Coulomb barrier and by 24% for the bombarding energy 1.05 times 
the value of the Coulomb barrier. Calculated complete fusion cross sections agree with the 
experimental data. The breakup probability of 9 Be decreases with decreasing bombarding 
energies towards the Coulomb barrier, with a hindrance factor of 0.66 — 0.76 for complete 
fusion which is similar to the experimental one of ~ 0.75 obtained for the studied reaction 
and to the experimental value of 0.68 ± 0.07 for the similar reaction 9 Be + 208 Pb. For a 
fixed bombarding energy, the breakup probability increases with decreasing orbital angular 
momentum and it becomes maximal for a central collision. 

The effect of the breakup of 9 Be on inelastic channel is similar to its effect on complete 
fusion channel. Inelastic cross sections decrease with decreasing bombarding energies. Re- 
sults obtained for E cm = 46 MeV qualitatively explain the experiment. The results are 
similar to those obtained when the double-folding potential determines the trajectories. 

In general, the results indicate that absorption channels for bombarding energies very 
close to the Coulomb barrier observed in experiment [|15[] can also be associated with orien- 
tation and static deformation effects of 9 Be which lower the value of the Coulomb barrier 
and increase the fusion cross section. Our results also support the molecular description of 
the highly asymmetric reaction 9 Be + 209 Bi in the framework of the two-center shell model. 
It would be interesting to study other reactions with 9 Be but a more general two-center shell 



model |23| would be needed in order to describe, e.g., arbitrary orientations of the intrinsic 
symmetry axes of deformed nuclei. Calculations for the studied reaction using a three-body 
model within a coupled discretised continuum channels (CDCC) formalism are in progress. 
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FIG. 1. Neutron levels for 9 Be + 209 Bi — ► 218 Fr as a function of the separation r between the 
nuclei. Levels are characterised by the total angular momentum projection j z on the internuclear 
axis: j z = 1/2 (dashed curves), j z = 3/2 (dotted curves), j z = 5/2 (dashed-dotted curves) and 
other values (solid curves). See text for further details. 
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FIG. 2. a) The angular momentum independent part of the total nucleus-nucleus potential as 
a function of the separation r between the nuclei. The adiabatic TCSM-potential with spherical 
nuclei and the sudden double-folding potential are shown by solid and dashed curves, respectively. 
The smooth macroscopic part of the adiabatic TCSM-potential is presented by a dotted curve. 
The adiabatic TCSM-potential for the prolate nuclei in the pole-to-pole configuration is shown by 
a dashed-dotted curve, b) The nuclear surface energy En with spherical nuclei is shown by a solid 
curve. The nuclear double-folding potential is presented by a dashed curve. See text for further 
details. 
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FIG. 3. Complete fusion cross sections df us as a function of the inverse of the bombarding 
energy 1/E c . m . in the center of mass system. Predictions (without breakup) of the Glas-Mosel 
model for the adiabatic TCSM-potential with spherical nuclei (see Fig. 2a) are shown by a solid 
curve and values taken into account the breakup of 9 Be, by a dotted curve. Experimental data 
for the studied reaction (full dots) are from and for 9 Be + 208 Pb (full triangles) from fi~2|] . 
Complete fusion cross sections aj us calculated with the Glas-Mosel model (without breakup) for 
the adiabatic TCSM-potential with prolate nuclei (see Fig. 2a) are presented by a dashed-dotted 
curve. See text for further details. 
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FIG. 4. Breakup probabilities Pbup as a function of the orbital angular momentum 1% (0 < I < ^ cr ) 
for captured trajectories obtained for the adiabatic TCSM-potential with spherical nuclei (see Fig. 
2a). Values for E c . m . = 42 MeV and E cm . = 46 MeV are shown by dashed and solid curves, 
respectively. See text for further details. 



21 




FIG. 5. Scattering function 9 c , m ,(l) (upper part), distance of closest approach r m - m of the nuclei 
(middle part) and dissipated energy E^iss (lower part) as a function of the orbital angular momen- 
tum Ifi (I > l cr ) for non-captured trajectories at E c . m . = 46 MeV. Values obtained for the adiabatic 
TCSM-potential with spherical nuclei (see Fig. 2a) and for the double-folding potential are shown 
by solid and dot-dashed curves, respectively. See text for further details. 
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FIG. 6. At E c _ m . = 46 MeV, dissipated energies Ediss (upper part) and inelastic classical cross 
sections dai n /dVt (middle part) as a function of the scattering angle 9 c .m. f° r inelastic non-captured 
trajectories obtained with the adiabatic TCSM-potential for spherical nuclei (see Fig. 2a) by 
using different values of the radial friction coefficient in units of 10 -23 MeV _1 s. For a radial 
friction coefficient c® = 4 x 10 -23 MeV -1 , the inelastic angular distribution taking into account the 
breakup of 9 Be is shown by a dashed-dotted curve. The experimental point (full dot) is from [15]. 
Elastic angular distributions da e i/daji (lower part), normalised to the Rutherford cross section 
ctr, obtained for the adiabatic TCSM-potential with spherical nuclei and for the double- folding 
potential are shown by solid and short-dashed curves, respectively. Experimental data (full dots) 



are from 1 151. See text for further details. 
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